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Decompositions of a Lie group G are methods to factorize every element leGas 

X = X 1 X 2 ---X m , (1.1) 

where the factors X±, . . . , X m belong to one dimensional subgroups. Decompositions of unitary evolutions in simpler 
• terms are of interest in quantum control and information theory for at least three reasons. They allow to simplify 
the task of controlling the evolution of a quantum system to a target into a sequence of simpler subtasks, consisting 
of control problems to more easily reachable targets (see e.g. 0,@>lp)- They allow the analysis of several features of 
quantum dynamics such as entanglement generation, time optimality and parameter identification (see e.g. 0, 0, 0- 
\Q , LlJj])- They give methods to produce unitary evolutions in a laboratory by combining a sequence of readily reproducible 
evolutions. In particular, in quantum information theory, a decomposition can be seen as a method to generate a 
quantum logic operation from a sequence of elementary operations ja, Q ■ 

In many cases, decompositions of the unitary group U (n) correspond to vector space decompositions of the cor- 
responding Lie algebra u(n) and each element H £ iu(n) represents the Hamiltonian of a possible evolution. In the 
analysis of multipartite quantum systems, it is useful to distinguish Hamiltonians acting on single subsystems, called 
local Hamiltonians, and Hamiltonians describing the coupling between two or more systems, called interaction ( or 
entangling) Hamiltonians. In the unitary group, these Hamiltonians generate local and entangling evolutions, respec- 
^ tively. In particular, if one considers a multipartite system composed of N subsystems of dimensions Ui,... , njv, 
^ , the space of all possible Hamiltonians is given by iu(niU2 • • • n^). Once we have an orthogonal basis in iu(rij), 
' j = 1, . . . , N, given by Hf. , lj = 1, . . . , nj, then a basis of u(nin2 ■ ■ ■ Un) is given by 

iHl®Hl®...®H? N . (1.2) 

The subalgebra of local Hamiltonians is spanned by elements where all the factors in the tensor product are equal 
to the identity except one. They produce unitary evolutions of the factorized form X\ g) X 2 <8> . . . (8) Xm which 
corresponds to transformations on the single subsystems. In this spirit, the decomposition given in [(| recursively 
factorizes a unitary transformation on n qubits in local and entangling transformations, the latter one acting on two 
subsystems at a time only. In a recent paper £| a method was given to generate decompositions in the tensor product 
space for general multipartite quantum systems of arbitrary dimensions starting from decompositions of evolutions 
on the single subsystems. 

In the spirit of the last two works cited, we present in this paper a recursive procedure to decompose the unitary 
evolution of a bipartite system of arbitrary dimensions so that every evolution is factorized into simple terms and it 
is clear what the entangling and local contributions of a single transformation are. The procedure we present applies 
recursively the Cartan decomposition of the Lie algebra so(n) by keeping the tensor product basis representation of 
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the Lie algebra u(n). In this basis, at the end of the procedure, it is easy to analyze the local and entangling character 
of each factor. 

The paper is organized as follows: in the next section we review some basic concepts and results on the Cartan 
classification of symmetric spaces and Cartan decompositions. We shall mention only the facts needed in the sequel of 
the paper. A detailed treatment can be found in [5|. The decomposition procedure we describe consists of an initial 
step which reduces the problem to a decomposition of the orthogonal group and a recursive procedure which allows to 
decompose further the elements of the orthogonal group. These steps are described in Sections llIII and lTVl respectively. 
We give several remarks highlighting the local and entangling factors needed in the factorization. Section[V]is devoted 
to a discussion and a numerical example. 

Notation: In the following we will use several times the definition of the elementary matrices E mn , 
and of their antisymmetric and symmetric superposition, respectively as 

^■mn • — E. mn E nrni ^l ran . — E mn -\- E nm . (1.4) 

We denote by A T the transposed of the matrix A. A n x m rectangular matrix is denoted by A nxm ; the n x n 
identity matrix by l n . Finally, we call a sign matrix a matrix of the form dia<?(±l, ±1, . . . , ±1) with all the possible 
combinations of + and — . 

II. BACKGROUND MATERIAL 

In the following, we shall use (in a recursive manner) Cartan decompositions of the Lie algebras su{n) and so{n) as 
well as Cartan decompositions of direct products of (isomorphic copies of) these Lie algebras. A Cartan decomposition 
of a semisimple Lie algebra [ is a vector space decomposition 

[ = «ep, (2.i) 

where 6 is a subalgebra, namely 

[e,t]ce, (2.2) 

and the complement subspace p = I x is such that 

[ P ,p]ce, [p,e]c P . (2.3) 

To a Cartan decomposition of I there corresponds a factorization of e , the connected Lie group associated to [, 
such that every element l£e' can be written as 

X = KP : (2.4) 

where K belongs to e { , the connected Lie group corresponding to 6, and P is the exponential of an element in p. The 
coset space e'/e { is called a symmetric space of e\ A maximal Abelian subalgebra of p is called a Cartan subalgebra 
associated to the decomposition and it is denoted by a. One can show, under appropriate assumptions, that 

|J KaK~ 1 =p, (2.5) 

Kee 1 

so that the factorization 12. 411 refines to 

X = K 1 AK 2 , (2.6) 

with Ki, K2 <E e { and A £ e a . The dimension of the Cartan subalgebra a is called the rank of the decomposition (or 
of the associated symmetric space). 

Cartan has classified all the symmetric spaces of the classical Lie groups, i.e. the Lie groups SU(n), Sp(n) and 
SO(n), and has shown that, up to conjugacy , the corresponding decompositions fall in one of few classes which he 
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has described. In particular, for [ = su(n), there are three types of decompositions labeled by AI, All, AIII. In the 
following, we shall use only decompositions of the type AI, given by 

t = so(n), p=so(n) ± . (2.7) 

The inner product in the space su(n) is given by (A, B) = Tr(AB*) so that so(n)- L is the subspace of su(n) spanned 
by purely imaginary matrices. The rank of this decomposition is n — 1. Decompositions of the type All and AIII will 
not be considered here. 

We consider now the Lie algebra I = so(n). If n = 1, this algebra contains only the null matrix 0. For n > 2, when 
n is odd there is only one type of Cartan decomposition, denoted by BDI. Fixing two positive integers r and q such 
that r 5* q 1 and r + q = n, the matrices k € t have the form 

k = (o °b) > ( 2 - 8 ) 

with A € so(r) and B £ so(q). Matrices p € p have the form 

P=(_ C T°), (2.9) 

for a general r x q matrix C . The rank of this decomposition is q. In this paper, we will not consider other 
decompositions. 



III. DECOMPOSITION OF UNITARY EVOLUTIONS IN C/(did 2 ); INITIAL STEP 



Consider two interacting quantum systems 5i and S 2 whose associated Hilbert spaces have dimensions d\ and d 2l 
respectively. According to the procedure described in Q, it is possible to obtain a decomposition for su{d\d 2 ) from 
decompositions of type AI of the Lie algebras associated to each subsystem, i.e. su(di) and su(c?2). We also include 
in the algebra scalar matrices, physically corresponding to shifts in the energy. We write 

u{dx) =so(di) eso(rfi)- 1 , u(d 2 ) = so{d 2 ) ®so(d 2 ) ± . (3.1) 

Let a J , j — 1,2, be a generic element of an orthogonal basis of iso(dj), and S"- 7 , j = 1,2, a generic element of an 
orthogonal basis of iso(dj)- 1 . Then the subalgebra of 11(^1^2) defined by 

t := spanker 1 ® S 2 , iS 1 <g> a 2 }, (3.2) 
along with its orthogonal complement in u(did2), 

p := spaniia 1 ® a 2 .iS 1 <g> S 2 }, (3.3) 
define a Cartan decomposition of u(did 2 ) as 

u(d 1 rf 2 )=t®p. (3.4) 

This decomposition is of type AI fj as t is conjugate to so(did 2 ) and p to so(did 2 ) ± . The rank of this decomposition 
is did 2 . A basis of the maximal abelian subalgebra o C p is given by tensor products of elements of the orthogonal 
basis of the maximal Abelian subalgebras associated to the single subsystems, which are of dimensions d\ and d 2 
respectively. Denoting by D 1 the diagonal elements of the type S 1 , and by D 2 those of the type S* 2 , a is given by 

a := span{iD x <8> D 2 }. (3.5) 

The associated Cartan factorization of A £ U{d\d 2 ) is 

X = KiAK 2 (3.6) 

according to the notation of the previous section. 
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Remark III.l Only in the simplest case of the decomposition of su(4) (i.e. d\ = 2 and d% = 2), studied for example 
in |10| . the decomposition 13. 6[) is a decomposition in local and nonlocal transformations. The local transformations 
are products of exponentials of matrices of the form iH ® 1 or 1 ® iH, where 1 is the identity matrix of appropriate 
dimensions and H is a generic matrix in ju(g?i) or i\x(d 2 ). Both local and nonlocal transformations are possibly 
present in the K\ and K 2 factors as well as in the A factor. However, obtaining a decomposition in terms of tensor 
product matrices will allow us to identify exactly where the local and nonlocal transformations are present in the final 
transformation. 

Remark III. 2 We notice that, in general, only one nonlocal transformation, along with the set of the local transfor- 
mations, is sufficient to obtain all the possible values for A in (|3.6|) . To see this, notice that the factor A is the finite 
product of exponentials of elements of the form ictjkEjj ®Ekk, with j = 1, . . . , dn k = 1, . . . , d 2 , and ctjk real numbers. 
Since, for every En is unitarily equivalent to En, the Hamiltonian H = En ®En, along with local transformations, 
is sufficient to generate any element of the form A. Notice that an alternative (universal) nonlocal Hamiltonian is 
given by an Ising interaction between two spins, which in our notation reads as (En — Ed^) <8> (En — Ed 2 d 2 )- 

We now turn our attention to decomposing the elements Ki and K 2 in 13. This will be obtained through a 
recursive procedure via iterate decompositions of so(n). 



IV. DECOMPOSITION OF UNITARY EVOLUTIONS IN U(d 1 d 2 ); RECURSIVE PROCEDURE 

The Lie algebra 6 defined in 13.211 is conjugate to so(d 1 d 2 ). We rewrite its definition below: 

t := spaniia 1 (g> S 2 , iS 1 <g> cr 2 }, 

with a° , j = 1,2, belonging to an orthogonal basis of iso(dj) and S J , j — 1,2, belonging to an orthogonal basis of 
iso(dj) ± . 

A special case arises when d\ = d 2 = 1, and only the matrix belongs to the corresponding Lie algebra. This case 
is not of physical interest as it would imply a one dimensional quantum system. However it may arise as the final step 
of the recursive procedure we are going to present. A special, nonphysical case is also di = 2 and d 2 — 1 or viceversa. 
In this case, the Lie algebra t contains only one element. Another special case is given by d\ = d 2 = 2. In this case, 
consider the Pauli matrices 



l\ _ (0 -A _ ( 1 

1 o ' ° v : ~ { i o r ° z : " U -1 



(4.1) 



and the 2x2 identity matrix 1. Then t = so(4) is the direct sum of two commuting subalgebras, Si and s 2 , each 
isomorphic to so(3), and given by 

51 :— span{ia y <E> l,icr x ® a y ,ia z ® a y }, 

5 2 := span{il ® a y ,io- y (g) o- x ,ia y (g) er z }. (4.2) 
Therefore K\ (and analogously K 2 ) in (|3.t)|l can be written as the product 

K l = F l F 2 = F 2 F 1 , (4.3) 

with F\ and F 2 in the Lie group corresponding to Si and S2 respectively. A Cartan decomposition can be performed 
on Si (and analogously on s 2 ) which is an Euler decomposition as Si is isomorphic to so(3), and allows to express F\ 
as 

Fi = LiNL 2 , (4.4) 

with Lj = e 03 " 7 "® 1 , j — 1, 2, for real parameters ctj, and N = e P l < , *® c 'y for a real parameter (3. Notice that L\ and L 2 
are local transformations while N is nonlocal. The same can be done for F 2 , moreover the nonlocal transformation 
for F 2 can be obtained using a local similarity transformation from the one for F\ or viceversa, so that for F\ and F 2 
we need only one nonlocal Hamiltonian. 

Consider now the case where at least one between d\ and d 2 is greater than 2. As it was done for the initial step in 
the previous section, we look for decompositions concerning the single subsystems to induce a decomposition on the 
total bipartite system. The elements of the type ia 1 and ia 1 are real, skew-symmetric, square matrices of dimensions 
d\ and d 2 respectively. Let d\ > 2 without loss of generality. On the Lie algebra of d\ x d\ skew-symmetric matrices 
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(that is, matrices of the type a 1 ) we perform a decomposition of the type BDI (see previous section) by selecting two 
positive integers r± ^ qi ^ 1 so that r\ + qi = d\. In the resulting Cartan decomposition, so(di) = t © p, the Lie 
algebra t is spanned by block diagonal skew-symmetric matrices with the upper block of dimension r\ and the lower 
block of dimension q\. We denote this type of matrices by ia 1,D (where D stands for 'diagonal'). The skew-symmetric 
matrices in the complement p will be denoted by ia 1,A (where A stands for 'anti-diagonal'). We as well separate 
matrices of the type S 1 into block diagonal and block anti-diagonal and denote them by iS 1,D and iS 1,A , respectively. 
Analogously, we define a decomposition of type BDI on so(c?2) introducing two positive integers r 2 ^ q 2 5^ 1, with 
r 2 + ?2 = d 2 and matrices of the type ia 2 ' D , ia 2 - A , iS 2 ' D and iS 2,A . In the special case where d 2 = 2, we can only 
choose r 2 = q 2 = 1 and we do not, in fact, obtain a decomposition of so{d 2 ) of the type BDI. However, we still 
formally decompose matrices of the form a 2 and S 2 in (block) diagonal and (block) anti-diagonal components and 
notice that, in this case, the only matrix of the type ia 2 ' D and iS 2 ' A is the 2x2 zero matrix. 

These decompositions on the two subsystems induce a decomposition on the overall bipartite system. More precisely 
we decompose 6 in l|3.2fl as follows 



t:= t'®p', (4.5) 



with 



f - spaniia 1 ' ® S 2 > D , iS 1 ' ® a 2 > D , ia^ A ® S 2 > A , iS^ A ® <r 2 > A }, 

p' = spaniia 1 ^ ® S 2 < D , zct 1 ^ ® S 2 " 4 , iS 1 ^ ® o 2 ' A , iS^ A ® a 2 ^}. (4.6) 

The following Theorem summarizes the features of this decomposition. It also gives, in its proof, a coordinate 
transformation to write the elements of the subalgebra 6' and its complement p' in the standard form. 

Theorem 1 The decomposition oft defined in ^~b\ ) is a Cartan decomposition, i.e. 

[t',f]ct', [e', P ']cp', [p', P ']cf. 

As a decomposition of so{d\d 2 ), it is a Cartan decomposition of type BDI with indices r and q satisfying r ^ q ^ 1, 
r + q = dxdz, and with 

r = nr 2 +qiq2, q = nq 2 +qir 2 . (4.7) 
Accordingly, the dimension of the associated Cartan subalgebra o'Cp' is q = r±q 2 + q\r 2 . 

Proof. We explicitly exhibit a conjugacy which transforms elements of t' into the form (|2.8I) and elements of p' into 
the form l|2.9|) . In particular, notice that the matrices ia 1,D ® S 2,D , iS 1,D ® a 2,D have the form 




\ 



(4.8) 



(4.9) 



-B 9lXgi <8> A, 2 x ?2 / 

while the matrices zcr 1,j4 ® S 2,A and iS' 1 '" 4 ® <7 2,j4 are of the form 

/ F riXgi ®G r2X(Z2 \ 

fr ._ ±F riXqi ®Gl Xr2 

T*£ X9l ®<W 

\--P£xr 1 ®G ! £xm 

A straightforward calculation shows that, defining 
the matrices 



nr 2 













































fci := Rk x R T , fc 2 := Rk 2 R T , 



(4.11) 
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have the form given in (|2.8|l where the upper block has dimension r = r\r 2 + 9192 and the lower block has dimension 
9 = r iQ2 + Qi r 2- Analogously one shows that the conjugacy defined in l|4.1U|l transforms elements in p' into elements 
of the form p in 1(2.9(1 ■ 

□ 

In view of the decomposition (|4.5|) any element K\ (and analogously for K 2 ) in 1(3.6(1 can be written as 

K x = K[A'K 2 , (4.12) 

where K[ and K' 2 belong to the Lie group associated to the Lie algebra 6', conjugate to so(r) © 50(9), and are to be 
further factorized. The matrix A 1 belongs to the Abelian Lie subgroup associated to the maximal Abelian subalgebra 
a' G p'. In the following Proposition we find an orthogonal basis for such a Cartan subalgebra expressing it in terms 
of tensor products 2 . 

Proposition IV. 1 The (ri +91)92 matrices 

Njk ■— Ejj © Ak t r 2 +k> j = !,■■■, n+Qi, k = l,...,q 2 , (4.13) 
along with the [r 2 — 92)91 matrices 

Mfr.= A Lri+ f®Eu, l = q 2 + l,...,r 2 , f = l,..., qi , (4.14) 
span a Cartan subalgebra a' 6 p'. 

Proof. The dimension of the vector space spanned by the matrices in (|4. 13|> , (|4.14(l is , in fact , (r 1 + 91 ) q 2 + (r 2 — q 2 ) 91 = 
^192 + 9i^2- Therefore, we only have to verify that matrices of the type 14.13(1 and 1(4.14(1 commute with each other. 
The commutator between two matrices of the type ((4.13(1 always vanishes. Analogously, matrices of the type ((4.14|) 
commute with each other. The Lie bracket of matrices of the type ((4.13(1 and 1(4.14(1 vanishes too, since the products 
of matrices A^ ^+fc and En in the second factors are always zero. □ 



Remark IV. 2 It follows from Proposition llV. ll that the element A' in 1(4.12(1 is the exponential of a linear combination 
of matrices ((4.13(1 and 1(4. 14|) or (equivalently) the product of exponentials of matrices proportional to these. The 
resulting unitary transformations may be entangling or local. However, since all the matrices of the form Ejj are 
unitarily equivalent to each other and the matrices of the type Aki are also unitarily equivalent to each other, only 
one (entangling) Hamiltonian of the type 1(4.13(1 . one of the type ((4.14|) along with local operations are sufficient (and 
necessary) to generate all the possible factors A' in 1(4.1211 . 

Remark IV. 3 A further reduction of the nonlocal Hamiltonians to be used is obtained by noticing that all the 
transformations in Remark lIII.2l and in Proposition IIV. II can be obtained with only one Ising Hamiltonian and local 
transformations. Therefore only one nonlocal Hamiltonian is needed to implement all of these transformations. 

We now further factorize the elements K [ and K' 2 in ((4.12(1 using, once again, a Cartan decomposition of the Lie 
algebra t' isomorphic to so(r) ©50(9), with r and 9 defined in ((4.7(1 . In particular, we decompose t' as follows 

t' = t"ffip", (4.15) 

with 

I" = span{i<j^ D © S 2 > D , S 1 ' 13 © i( r 2 ' D }, 

p" = span{ia^ A © S 2 ' A , S^ A © ia 2 < A }. (4.16) 
The matrices in {" are block diagonal matrices and t" is 

t" = Bo{r x r 2 ) © so(ri9 2 ) © 50(91 r 2 ) © 50(9192), (4.17) 



2 An alternative procedure is to transform the Lie algebra 6' according to the change of coordinates 14.101 and to find the Cartan 
subalgebra in the standard basis 
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where each term refers to a block on the diagonal. For example, the first block, corresponding to So(rir2), contains 
matrices obtained as tensor products ia 1,D <g> S 2 ' D or S ' ® cr 2 ' D where all matrices involved have the second block 
equal to zero. This corresponds to two decompositions of the type BDI, one on so(r) and the other on so(q). The 
Cartan subalgebra in p" is the direct sum of the two Cartan subalgebras of the two decompositions. It has dimension 
qiQ2 + min{riq2, q\r2} ■ The following proposition explains how to find a basis of this Cartan subalgebra as tensor 
product matrices 3 . 



Proposition IV. 4 A Cartan subalgebra of the decomposition (|4.15(l . 14.16fl is spanned by the (71(72 matrices 



JV, 



E 



j 1 



^ m ,in--y_ E T j®y E i n 

with 1 < J < r*i, 1 ^ I ^ q% and 1 ^ m ^ ri , 1 ^ n ^ (72 satisfying 

(j - 1)^2 + m = s, (I - l)q 2 + n 
with s = 1, . . . , (71(72, along with matrices 



E n 



Eji 
El 



E„ 



-El 



0. 



M 



jm,ln ■ — 



En , 
-ET t I 



E m , 

f t n 

ran 



E jt 
El 



(81 



E, 



-El 



run 





with 1 ^ j ^ r\, 1 / ^ q\ and 1 $5 m ^ r2, 1 ^ n ^ Q2 satisfying 

(j - 1)92 + n = s, (i - l)r 2 + to = s 

with s = 1, . . . , min{riq2,qir2}- 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



Remark IV. 5 For all s, there corresponds a unique pair (l,n) so that the second relation in (|4.19|) is verified. There 
is some freedom in choosing the pairs (j, to) satisfying the first relation in 14.19f) . However, for every value of s, and 
therefore of / and n, one is allowed to choose a unique pair (j, to). An analogous meaning has the notation in (|4.2()(l 
and |g2D>. 

Proof. Matrices of the form (|4.18(l commute with matrices of the form (|4.2()JI . since these matrices form Cartan 
subalgebras associated to decompositions of so(r) and so (q), respectively. To show that matrices of type (|4.18|) 
commute, one verifies that the commutators of two matrices corresponding to indices (jimi,lini) and (j2"^2, ^2^2) 
vanish. In fact all the blocks of such matrices are zero except for the 1,1 and 2,2 blocks which, from a direct 
calculation, turn out to be equal to 



^{Ej 2 jiSi 1 i 2 E m2rni 5 nin2 Ej 1 j 2 S[ 1 i 2 E miJn2 5 nin2 ), 



and 



<E l2h 6 jlh ®E. 



(4.22) 



(4.23) 



respectively. However these are also zero if n\ ^ «2 and/or l\ ^ l 2 as well as in the case 1% = I2, rii — n 2 (and 
therefore j\ = j 2 , mi = m 2 , see Remark IIV.5|) . A perfectly analogous argument holds in the case of commutators of 
matrices of the form (|4.20() . □ 



Remark IV. 6 Notice that all the Hamiltonians l|4.18|l are locally unitarily equivalent to each other. The same is 
true for the Hamiltonians (|4.20l) . Therefore only two more entangling Hamiltonians are needed. 

At this point we are left with the Lie algebra so (pipi) © 30(51(72) ©50(^1(72) ffiso(p29i)- The construction proceeds 
recursively by decomposing each one of the four component Lie algebras and so on, until one finds one of the Lie 
algebras so(l) (which we define as the element zero), so(2) (which consists of a single element), so(3) or so(4) (which 
are treated as it was explained at the beginning of the procedure) . 



Alternatively, one can construct a basis for this Cartan subalgebra working in the standard representation using the coniugacy given in 
fTTul 
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V. DISCUSSION AND A NUMERICAL EXAMPLE 



In order to illustrate the Lie group decomposition described in the previous sections, we consider the generalized 
SWAP operator X sw acting on three qubits and rotating their states in a cyclic fashion. Its action is defined in the 
tensor product basis as 



|*)i O \j) 2 <8> \k) 3 -> <8> \ih ® |i) 3 , 



(5.1) 



where i, j, k — 0,1 and {|0}, l)}i,2,3 are orthonormal basis for the Hilbert spaces of the three systems. 

This operator is relevant in quantum information and computation since it enables to switch the quantum states 
of different systems. For example, assume that one is interested in the state of the third system, but only the first 
system is accessible and can be controlled; then the application of the generalized SWAP operator will enable to 
transfer the state of the third system to the first system. As local operations alone clearly cannot implement the 
generalized SWAP and this has to involve some degree of entanglement among the various subsystems. We consider 
an hypothetical situation where it is possible to create interaction between the first qubit and the other two as a 
whole although it is difficult to create interactions with the single qubits 2 and 3. This justifies to consider the total 
Hilbert space as the tensor product of a 2-dimensional subspace with a 4-dimensional one (that is, d\ = 2 and d 2 = 4). 
Therefore X sw 6 U(8) will be decomposed accordingly 4 . 

In the specified basis, with standard ordering, the matrix representation of this operator is given by 

/l 0\ 
1 
1 
1 
1 
1 
1 
\0 1/ 

This transformation belongs to 50(8), therefore the first step of the decomposition of U(8) is trivial: Ki = X sw and 
A = K 2 = 1. 

For the first step of the recursive part of the procedure, we choose ri = qi = 1 and r 2 = q 2 = 2. We find it 
convenient to work in the basis of the Hilbert space such that t' — so(2) © so(4), obtained by performing the change 
of basis R\ijk) — > \ijk)' , with R given in l|4.10|l . which, in this particular case (ri = q\ = 1, r 2 = q 2 = 2), takes the 
form 



(5.2) 



R = 



'1 2 0' 

1 2 

1 2 

1 2 0. 



(5.3) 



In these coordinates, the SWAP operator is written as X sw — RX SW R T , that is 



A 


































1 
































1 











1 
































1 











1 




















1 




















Vo 














1 





0/ 



(5.4) 



The elements of the Cartan subalgebra a' (defined in Proposition II V . 1 \ . which in this case are only of the form 14.13fl . 
are transformed by R into elements of the form 

' D 
-D 



4 We believe that extensions of the procedure presented here to multipartite systems are possible at the price of an increased notational 
complexity. 
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The computational problem is to find two 4x4 diagonal matrices, D x and D 2 , with 

A' = e 5 ' 

A' e SO(8), and matrices £ SO (4), i,j = 1,2, such that 



D 1 D 2 
~D 2 D 1 



1 = I Ky 



K' 2 = 



K21 
K z 



(5.5) 



and = K' 1 A'K 2 . To perform this task we propose an algorithm which uses ideas similar to the ones for other 
Cartan decompositions (cf. e.g. and the references therein). We illustrate this algorithm for the dimensions of our 
problem but generalizations to other dimensions are obvious. Let us write X sw with 4x4 blocks X«, i, j — 1, 2, as 



X sm — 



Xn Xi 2 

X 2 i x 22 . 



Equation l|5.5|l is equivalent to the four matrix equations 



X\ 2 

X 2 \ 

x 22 



K XX D X K 2X , 
K 1X D 2 K 22 , 
-K X2 D 2 K 2X , 
K 12 D X K 22 . 



From the first one, we obtain 



XwX^Kyy — K XX D\, 



(5.6) 



(5.7) 
(5.8) 
(5.9) 
(5.10) 



(5.11) 



which is an eigenvalue equation as D\ is diagonal. In the generic case, when all the eigenvalues oi XnX^ are different, 
equation (|5.11|l determines Ku and D\ up to the right product by a sign matrix (and the fact that det(Kn) = 1). 
Moreover, D x gives D 2 up to a sign matrix from the requirement that D\ + D\ = I4. Using K X i and D 2 in l|5.8|l . we 
obtain K 22 up to a sign matrix. Plugging K 22 in l|5.10|l . we get K X2 up to a sign matrix and from (|5.9() we find ^21- 
Finally, we adjust the sign matrices to make l|5.7|l through (|5.1U|I consistently verified. In the case where X^X^ has 
multiple eigenvalues, there is more freedom in the choice of K xx at the initial step, but then one proceeds in the same 
way and determines the other matrices up to some degree of freedom. At the end of the procedure, these degrees of 
freedom are exploited to make (|5.7|I - H5.10|I jointly satisfied. 



Using this procedure, we have found for our example 



K xx = 



K 2X 



/l 
10 
10 

Vo -u 

/l 0\ 



K X2 =- 



1 





1 



D x 





Vo 1 

(\ ( 
( 
( 

Vo - 



Kn 



/0 IN 
10 
10 

Vi 0/ 

-1 0\ 

1 



1 0/ 

/o 0\ 
10 
10 

0000/ 



(5.12) 



We can repeat the same procedure as above to further factorize K xx , K X2 , K 2X , K 22 , and to obtain K[ = K"A'{K 2 
and K 2 = K'^A'^K". We finally get the decomposition of X sw , 



X sw = K'{A'[K'iA'K'iA'lKl 



(5.13) 
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where 



A'{ = diag{ 



A>> 



Vo 

diag{l4,, 



(\ 0\ 

1 

10 

1 0; 



'1 0> 

1 

10 

v0 -1 0/ 



fQ -1 0\ 

1 

10 

vO -1 0/ 



}• 



K'{ 



'0 -1 
.1 o 



-1 

1 



1 

-1 



}• 



K" = 



'0 -1 
,1 



l2,l 2 }, 

-1 N 





1 



dmg{l 2 , 
diag{l 2 , 
diag{l 2l I2, 12, 
diag{l 2 , 12, 

We now use the transformation R in l|5.3|l to write X sw in the original coordinates as 

X sw = K'lAlK'iA'K'lA'iK'l, 



(5.14) 



-1 

1 



}• 

,1 2 }- 



(5.15) 



with A' = R T A'R, A" = R T A"R and K'j = R T K'-R 1 j = 1, 2. We can write all the factors in (|5~T5|) as exponentials 
of appropriate matrices in the tensor product basis: 

A' = e a 



K" 

j 



e a ", fc = l,2, 



(5.16) 



where 



a = 




A 24 ) - 


37T 

f y(#22 <8 A13) +7T(^22 <8 


A24), 
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a i = 


-(A 12 € 


) ^24 + 


Ql2 <8 A 24 ) + ^(Al2 <8 ^24 


- O12 (8) A 2 4 


n 

a 2 = 


^(Ai2 <8^24 - 


^12 (8 A 24 ), 




k'( = 


"^-(^22 


8> A 34 ) 


7T . . 37T , 

+ -(^22 C8 A 12 ) + -yC®" 


8 A 34 ), 


K = 


■y (-^22 


8 A 34 ) 


1 




fc 3 = 


■y (-^22 


8 A ia ) 


) 




&4 — 




8 A 34 ) 







37T 

— (Ai2<8)fii 3 -fiia® Ai 3 ), 



(5.17) 



It is interesting to observe what number of nonlocal transformations are needed to perform the given task if we 
are able to perform any local transformation on the two subsystems. Notice that we are considering the system 
as a bipartite system of a two level system with a four level system. In essence, we assume that we have to decide 
appropriate interactions between the two subsystems (two an four dimensional) which along with local transformations 
will allow us to perform the given task. By grouping the matrices that are equivalent through local similarity 
transformations it is clear that the Hamiltonians 



H x := £n®A 3 4, 

H 2 := A12 (8 ^24 + fiia <8 A 24 , 

H 3 := A12 ® ^24 - H 12 ® A 24 



(5.18) 
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are sufficient 5 . 
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5 In view of 15.171 one can replace H2 with H2 '■= A12 ® Q.24, 



